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ABSTRACT 

We present an analysis of the mass distribution in the core of A 383 (z = 0.188), one of twelve X-ray luminous 
galaxy clusters at z ~ 0.2 selected for a comprehensive and unbiased study of the mass distribution in massive 
galaxy clusters. Deep optical imaging performed by the Hubble Space Telescope (HST) reveals a wide variety of 
gravitationally lensed features in the core of A 383, including a giant arc formed from the strongly-lensed images 
of two background galaxies, two radial arcs in the halo of the central cluster galaxy, several multiply-imaged 
arcs and numerous arclets. Based upon the constraints from the various lensed features, as well as from color 
information from ground-based observations, we construct a detailed model of the mass distribution in the central 
regions of the cluster, taking into account both the cluster-scale potential and perturbations from individual cluster 
galaxies. Keck spectroscopy of one component of the giant arc identifies it as an image of a star-forming galaxy at 
z = 1.01 and provides an accurate measurement of the mass of the cluster within the projected radius of the giant 
arc (65 kpc) of (3.5 ±0.1) x 10'^ M©. Using the weak shear measured from our observations we extend our 
mass model to larger scales and determine a mass of (1.8±0.2) X 10'"* M© within a radius of 250kpc. On smaller 
scales we employ the radial arcs as probes of the shape of the mass distribution in the cluster core (r < 20kpc), 
find that the mass profile is more peaked than a single Navarro, Frenk & White (NEW, 1997) profile. Our findings 
therefore support the proposal that massive cluster cores contain more mass than can be explained by a single 
cluster-scale NEW profile. The optical and X-ray properties of A 383 indicate the presence of a central cooling 
flow, for which we derive a mass deposition rate of > 200 M© yr"' . We also use the X-ray emission from A 383 to 
obtain independent estimates of the total mass within projected radii of 65 and 250kpc: (4.0^[ ') x 10'^ M© and 
(1 .2 ± 0.5) X lO'^^M©, which are consistent with the lensing measurements. 

Subject headings: cosmology: observations — gravitational lensing — clusters of galaxies: individual: A 383 — 
galaxies: evolution 



1. INTRODUCTION 

Massive clusters of galaxies represent one extreme of the 
mass spectrum of collapsed structures at the present day. Their 
properties are expected to reflect predominantly gravitational 
processes and therefore provide unique insights into the nature 
and distribution of the dark matter which drives the formation 
of structure (e.g. Eke et al. 1996; Viana & Liddle 1996; Bahcall, 
Ean & Cen 1997; Kay & Bower 1999). In particular, accurate 
measurements of the mass distribution in clusters across a range 
of scales can be used to test the claim of a universal form for 
the profiles of dark matter halos (Navarro, Erenk & White 1997, 
NEW) and hence the nature of dark matter. At large scales the 
NEW profile falls off as p oc r"^, steeper than an isothermal 
model (p (X r~^). On smaller scales the NEW profile breaks to 
a shallower slope with p oc whilst retaining a central cusp. 
Analyses of the mass profiles of dwarf galaxies from dynami- 
cal studies have used the apparent lack of a central cusp in the 
density distribution to reject the form of the NEW profile and 
argue instead for the existence of self-interacting dark matter 
(e.g. Moore et al. 1998). However, on larger scales, (e.g. lu- 
minous elliptical galaxies and massive clusters of galaxies) the 
form of the mass profile has yet to be investigated systemati- 



cally. 

The central regions of massive, compact clusters at moder- 
ate redshifts (z < 0.5) act as strong gravitational lenses forming 
multiple images of serendipitously placed background galaxies 
(e.g. Small et al. 1996). The properties of the images of these 
background galaxies (position, relative surface brightness and 
parity) can be used to accurately model the distribution of to- 
tal mass (both baryonic and non-baryonic) within the cluster 
core (Kneib et al. 1996, K96; Small et al. 1996; Natarajan et 
al. 1998). Such mass maps provide the most direct and detailed 
view of the distribution and morphology of dark matter in the 
centers of galaxy clusters. Strong gravitational lensing by rich 
clusters is therefore an ideal observational tool with which to 
test the halo properties predicted by cosmological models. 

The first mass maps constructed from lensing observations 
were produced in the early and mid-1990's (Tyson et al. 1990; 
Meflier et al. 1993; Eahlman et al. 1994; Kneib et al. 1994; 
Small et al. 1995). A significant improvement arrived with the 
refurbishment of HST which provided the resolution necessary 
to identify faint, lensed features in the crowded cores of rich 
clusters (e.g. K96), thereby allowing more detailed mass mod- 
els to be constructed. 
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On small scales, unique constraints can be obtained from ra- 
dial arcs (rare radially magnified images of background galax- 
ies) which are found in the very centers of a few cluster lenses. 
Such features are very difficult to detect in ground-based ob- 
servations, with only a single example known (MS 2137.3-23; 
Fort et al. 1992; MelHer et al. 1993). In contrast to this, HST has 
uncovered several radial arcs in previously well-studied cluster 
lenses (e.g. A 370, Small et al. 1996; AC 114, Natarajan et al. 
1998) which enabled the first tests of the form of cluster mass 
profiles on scales of less than ^ 100 kpc (Williams et al. 1999). 

Previous HST studies of rich clusters covered a heteroge- 
neous mix of clusters, either selected because they were pre- 
viously well-studied (e.g. Small et al. 1997) or because they 
were known to be strong lenses (e.g. K96; Small et al. 1996). 
To study the form of the mass profile in rich clusters in an un- 
biased fashion we need i/^Jobservations of an objectively se- 
lected cluster sample. Ideally such a sample would be mass- 
selected; however, in the absence of samples compiled from 
large-scale weak lensing surveys. X-ray selected cluster sam- 
ples (Gioia et al. 1990; Ebehng et al. 1998, 2000; DeGrandi 
et al. 1999; Ebeling, Edge & Henry 2000) are best suited to 
selecting well-defined samples of massive clusters. 

Following this premise, we are conducting a survey of twelve 
of the most X-ray luminous clusters (Lx > 8 x 10"*^ erg s~' , 0. 1- 
2.4 keV) in a narrow redshift slice at z ^ 0.2, selected from the 
XBACs sample (X-ray Brightest Abell-type Clusters; Ebehng 
et al. 1996). As XBACS is restricted to Abell clusters (Abell, 
Corwin & Olowin 1989), it is X-ray flux limited and not truly 
X-ray selected. However, a comparison with the X-ray selected 
ROSAT Brightest Cluster Sample (BCS, Ebehng et al. 1998, 
2000) shows that - 75% of the BCS clusters in the redshift and 
X-ray luminosity range of our sample are in fact Abell clusters. 
Hence, our XBACs sample is, in all practical aspects, indistin- 
guishable from an X-ray selected sample. 

In this paper we describe HST observations of the core of one 
of the first clusters observed in our survey, A 383 (z = 0.188), 
which reveal a multitude of strongly-lensed features. In com- 
bination with color information obtained from ground-based 
imaging observations in three passbands, the lensed features 
are used to produce an accurate model of the mass distribution 
within the central 500 kpc of the cluster We complement these 
results with an analysis of the X-ray properties of A 383 using 
archival ROSAT HRI data. In §2 we describe the observational 
data, their reduction and analysis, followed in §3 by a discus- 
sion of the interpretation of these data and the construction of 
the lens model of the cluster Finally in §4 we summarize the 
main results of our analysis and present our conclusions. We 
adopt Hq = 50kms"' Mpc"^ and qo = 0.5 throughout. At the 
cluster redshift 1" = 4.0 kpc in this cosmology. 

2. DATA REDUCTION AND ANALYSIS 

A 383 is a massive cluster of galaxies at z = 0. 188 with a core 
dominated by a luminous cD galaxy (Figs. 1 & 3). The X-ray 
luminosity from a pointed observation with the ROSAT HRI is 
(9.8±0.3) X lO^^ergs-i in the 0.1-2.4keV band (see §2.4). 

2.1. HST Imaging 

A 383 was observed wifli the HST WFPC2 camera on Jan- 
uary 25, 2000. Three exposures totaling 7.5 ks were taken 
through the F702W filter. Each exposure was shifted rela- 
tive to the others by 10 WFC pixels (1.0") providing a partial 
overlap of the chip fields. After pipeline processing, standard 
IRAF/STSDAS routines were employed to shift and combine the 



frames to remove both cosmic rays and hot pixels. Corrections 
for under-sampling of the point spread function and geometric 
distortion of the optics were made using the dither package 
within IRAF (Fruchter & Hook 1997). The final frame (Fig. 1) 
has a pixel scale of 0.05", an effective resolution of 0.15" and 
aid detection limit within the seeing disk of Rj()2 ~ 3 1 . 

To produce a catalog of faint arclets from our data we first 
analysed the HST frame using the SExtractor package (Bertin & 
Amouts 1996). All objects with isophotal areas in excess of 30 
pixels (0.074 arcsec ) at the /i702 — 24.5 mag arcsec isophote 
(3cr/pixel) were selected. A steep roll-over in the observed dif- 
ferential number counts arising from incompleteness occurs at 
R-JQ2 ~ 25.5. We estimate the completeness to be ~ 80% (ap- 
proximately 5(t) at Rjo2 ^ 25.5, which we adopt as our mag- 
nitude limit, giving a total of 457 sources within a 4.9 arcmin^ 
area (excluding the PC chip). Seventeen of these sources are 
classified as star-like on the basis of their profile shapes and are 
excluded from our analysis. 

For the purposes of the photometric analysis presented in Ta- 
ble 1 we correct the F702W photometry (^^702) to the Cousins 
/J-band in the following manner. The observed range in (B-R) 
colors of the arcs roughly translates to 0.2 < (V-R) < 1.0 giv- 
ing a typical color of (V-R) ^ 0.6. Adopting this color we 
obtain a correction of R-1Q2-R = 0.2 (Holtzman et al. 1995), and 
estimate that this correction introduces a systematic uncertainty 
in the /?-band photometry of ±0.06 magnitudes. 

The final HST frame (Fig. 1) reveals many previously un- 
known strongly lensed features, including a giant arc and two 
radial arcs, making A 383 a striking new addition to the catalog 
of cluster lenses at intermediate redshifts. 

2.2. Ground-based Imaging 

With the aim of extending our lensing mass map to the turn- 
around radius of the cluster, we have used the 3.6m Canada- 
France-Hawaii Telescope (CFHT) with the CFH12K camera to 
obtain panoramic images of A 383 on the nights of November 
14-16, 1999. Total exposure times of 7.2ks, 6.0ks and 3.6 ks, 
accumulated at 6-10 dither positions, were acquired in the fi, /? 
and / bands respectively. Data reduction was performed within 
IRAF using the MSCRED package including standard bias sub- 
traction and flat-fielding using twilight flats. The dithered ex- 
posures were aligned with the Digital Sky Survey frame of the 
same field to an rms accuracy of 0.15". More information on 
the reduction and analysis of these observations will be pre- 
sented in a forthcoming paper (Czoske et al. 2000). Here we use 
the central regions of the B and / band frames to provide pho- 
tometry of the lensed features in the cluster core (see Table 1). 
These two frames have seeing of 0.88" and 0.71" FWHM re- 
spectively. 

In addition to the CFH12K imaging, we have obtained K- 
band images of the core of A 383 with the UFTI imager on 
the 3.8m United Kingdom Infrared Telescope (UKIRT), Mauna 
Kea, on October 14, 1999. The final frame was accumulated in 
27 dithered sub-exposures of 90 s duration each to give a total 
on-source integration time of 2.4 ks, all in photometric condi- 
tions. Employing standard procedures these frames were re- 
duced, combined and calibrated using observations of UKIRT 
faint standards bracketing the science exposures. The final 
frame has seeing of 0.42" FWHM with 0.09"/pixel sampling 
and an effective 5a depth ofK = 20.3. 

Optical and optical-infrared colors of the lensed features in 
A 383 were measured off the ahgned and seeing-matchedB/?/^r 
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Fig. 1 . - The central region of A 383 as observed with HST WFPC2, overlaid with an isodensity representation of our lens model. 
Contours correspond to projected surface mass densities of 3, 4, 6, 8, 11, 15 x IO^Mq kpc~^. The numerical labels indicate the 
cluster members used in the lens model. The alpha-numeric labels identify the multiple images used to constrain the lens model and 
a number of other singly-imaged arclets. The spiral galaxy labeled B16 was included in our spectroscopic sample and has z = 0.6560. 
The disturbed morphology of the central galaxy, including what appears to be a dust lane passing between two peaks in the Ught 
profile, is apparent in the center of this figure. 
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frames using apertures designed for individual arcs. The arcs 
were first masked out of the science frame and the sky back- 
ground estimated by median smoothing over the aperture de- 
fined by the absent arc. The sky subtraction was then performed 
and the individual arc apertures applied to the sky-subtracted 
frame to obtain the BRIK photometry presented in Table 1 . Er- 
ror bars were estimated by varying the smoothing length used 
in estimating the sky background. 

2.3. Spectroscopy 

On January 26, 2000 we observed the cD galaxy in A 383 
with the Wide Field Grism Spectrograph on the University of 
Hawaii's 2.2-m telescope. The low dispersion spectrum ob- 
tained in a 1.8 ks exposure in extremely poor seeing (2.7") 
shows the full range of emission lines typically found in cD 
galaxies in cooUng flow clusters (e.g. Crawford et al. 1999) and 
places file cD at a heliocentric redshift of z = 0. 1880 ± 0.0012. 

On January 29, 2000 we obtained 28 additional spectra with 
the LRIS spectrograph (Oke et al. 1995) on the Keck-II 10-m 
telescope in MOS mode, with a total exposure time of 3.6 ks in 
average seeing of 0.9". Use of the 300/5000 grating centered 
at 7500 A provided wide spectral coverage (5000-10,000 A) 
at a spectral resolution of 2.55 A/pixel. Redshifts were mea- 
sured independently by three members of our team (OC, HE, 
JPK) and also using the RVSAO cross-correlation package un- 
der IRAF to accurately estimate the errors. Successive 3(t clip- 
ping around the peak of the radial velocity distribution yields 
a redshift of z = 0. 1880 ± 0.0002 for A 383 and a velocity dis- 
persion of cr = 1150;!;2oo kms"' from 18 galaxy redshifts. Ten 
other galaxies were found to be background, including B16 at 
z = 0.656, which lies close to the cluster core (see Fig. 1). 
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Fig. 2. - The spectrum of arc BOa in A 383 taken with the LRIS spectrograph on the 
Keck-II 10m telescope. The arc exhibits a blue continuimi and we identify the strong emis- 
sion line as [On]A3727 and confirm this with several UV absorption features (marked) to 
give a redshift of z = 1.0103 ± 0.0001. The upper spectrum is smoothed to tlie nominal 
resolution of the spectrograph and offset vertically for clarity. The shaded regions show 
areas of the spectrum strongly affected by night sky lines. 

We also obtained a spectrum of the lensed feature BOa which 
we show in Fig. 2. We identify the strong emission line at 
7493A as [Oil] A3727 which places the galaxy at a redshift 
of z = 1.0103 ± 0.0001. This interpretation is confirmed by 
the identification of Fell A2600 and Mgll A2800 absorption fea- 
tures. The properties of this galaxy are discussed further in §3. 

2.4. ROSAT/HRI Observation 

A 383 was observed with the ROSATHRI in February 1995 
for a total integration time of 27.4 ks (ROR #800890). An adap- 
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lively smoothed map of the observed X-ray emission (created 
using the ASMOOTH algorithm of Ebeling, White & Rangara- 
jan 2000) is shown in Fig. 3. We used a modified version of 
Steve Snowden's CAST hri software to correct for exposure 
time variations across the field of view and subtract a particle 
background component. 

From the fully processed ROSAT HRI image the total 
background-corrected///?/ count rate from A 383 is measured 
to be (0. 128 ± 0.005) ct s'^ within 1 .5 Mpc (6.2') of flie cluster 
center, using the mean diffuse background at r > 2 Mpc (8.2')- 
Assuming a standard plasma emission spectrum, a Galactic hy- 
drogen column density of 4.1 x 10^" cm"^ for this sight-line 
(Dickey & Lockman 1990), a metallicity of 0.3 and an ambient 
X-ray temperature of ~ 7.1 keV (estimated from the cooling- 
flow-corrected cluster Lx-kT relation of Allen & Fabian 1998) 
we derive a total unabsorbed cluster flux of (6.55 ± 0.23) x 
10"^^ ergcm"^ s~^ (0. 1-2.4 keV). The corresponding X-ray lu- 
minosity is (9.8 ± 0.3) X 10"^ erg s"' (0.1-2.4 keV), in agree- 
ment with the XBACs value of (8.0 ± 2.4) x lO'^ erg s'^ in the 
same energy band (Ebeling et al. 1996). 



fig3.gif 



Fig. 3. - X-ray flux contours of the adaptively smoothed emission from A 383 as seen 
with the ROSAT HRI. overlaid on an /?-band image of the cluster taken with the CF11I2K 
camera (see §2.2). All X-ray features are at least 3(t significant with respect to the local 
background. The lowest contour lies 30% above the background, adjacent contours differ 
by 30%. We adjusted the astrometry of the HRI image by 4" (well within the errors of the 
astrometry solution) to align the X-ray peak with the cD galaxy. The near circular symme- 
try of the X-ray flux contours reflects the symmetry of the lensing mass contoixrs plotted 
in Fig. 1. The region shown in our Fig. 1 is marked by the dashed box. 

3. MODELLING AND RESULTS 

In this section we first describe and interpret the multiply- 
imaged features observed in our HST frame of A 383 and ex- 
plain how these were used to constrain the mass distribution in 
the central regions of the cluster. We then discuss the form of 
the mass profile we derive in the very center of the cluster and 
compare this to theoretical expectations from high resolution 
N-body simulations. Finally we discuss the X-ray properties of 
A 383 and compare the mass distribution determined from the 
lens model with the estimates based on the ROSAT diA&. 

3.1. Multiply Imaged Features 

As a result of the superb resolution of our HST images, a 
detailed analysis of the morphology of the observed arcs al- 
lows us to identify seven multiply-imaged systems in the field 
of A 383. This is significantly greater than the number typically 
seen in ground-based observations of clusters. We discuss each 
multiply-imaged system below, and summarize their photomet- 
ric and spectroscopic properties in Table 1. The lensed features 
are identified in Fig. 1 and the multiple images are further illus- 
trated below in Fig. 4-6. 
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3.1.1. The Giant and Radial Arcs 

We interpret the giant and radial arcs in A 383 as being 
the lensed images of two background galaxies, BO and Bl, at 
z= 1.01 (§2.3) and z ^ 1.1 respectively. This interpretation 
arises from modelling (§3.2) of the radial arcs, which constrains 
the outer radial arc (Bid) to be a counter image of Bla/b/c at 
z ~ 1.1, and the inner radial arc (BOb) to be a counter image of 
BOa. BOa is also slightly closer to the center of the cluster than 
other components of the giant arc, supporting the suggestion 
that it is at a slightly lower redshift than B 1 . 




Fig. 4. - The giant arc in A 383 (see text for interpretation). Eacli tick mark represents 
1", and the orientation is as in Fig. 1. The dotted lines indicate the path of the z = 1.1 
critical line close to the observed arcs. 

BO (Figs. 4 &5) — This system comprises the brightest com- 
ponent of the giant arc (BOa) and the inner of the two radial 
arcs (BOb). It is also plausible that B15 is a component of the 
BO system given its proximity to BOa in the image plane, how- 
ever it's redder color (Table 1) suggests that it is either singly 
imaged, or that BOb has a significant color gradient. We are 
unable to measure the color of BOb to confirm this possibility 
due to the light from the central galaxy. We have also estimated 
the unlensed R-hand magnitude of BO (Table 1) which, along 
with it's blue optical-infrared color, suggests that BO is a low 
luminosity star-forming galaxy (My ~ -19.3). 




Fig. 5. - a view of the radial arcs as seen on the HST frame (enhanced by subtracting a 
median-smoothed frame from the science frame). Each tick mark represents I", and the 
orientation is as in Fig. 1. 

Bl (Figs. 4 & 5) — Bla/bis a pair of merging images strad- 
dling the z ~ 1.1 critical hne with counter-images Blc, at the 
extreme right hand end of the giant arc and B 1 d, the outer of the 



two radial arcs. The lens model (§3.2) and photometric analy- 
sis (Table 1) support this interpretation and appear to exclude 
B15 from being a further counter-image. The unlensed /?-band 
magnitude of Bl has been estimated (Table 1), from which we 
determine that it is a very low luminosity galaxy (My ~ -18.3). 

3.1.2. Perturbations by Galaxy Halos 

B2 andB3 (Fig. 6) — The lensing effect of neighboring clus- 
ter galaxy halos on this complex group of blue arcs demon- 
strates how galaxy-scale masses measurably perturb the cluster 
potential. We have modeled the effect of these perturbations 
and resolved B2 and B3 into two background galaxies, both at 
a redshift of z ~ 3.5. Although the photometric analysis is ham- 
pered by the nearby galaxy halos, the data presented in Table 1 
lend additional support to this interpretation. 

B2 comprises five images, of which B2a/b and B2d/e are 
pairs of images straddling the z 3.5 critical hne as it passes 
around the halos of cluster ellipticals #46 and #23 respectively. 
B2c lies outside this critical line and has the same parity as B2a. 
A sixth candidate counter image lies approximately 1" to the 
left of B3b under the halo of cluster elliptical #23. Multi-band 
HST imaging is necessary to achieve the accurate photometry 
required to confirm this hypothesis. There is also a low signal- 
to-noise feature lying in the saddle region between cluster el- 
lipticals #46 and #76, which is probably an extension of B2c. 
Counter images of this portion of the B2 system are undetected, 
due to the intrinsic faintness of the galaxy and the relatively low 
magnification of the counter images relative to the influence of 
the saddle on the observed feature. Finally, we estimate that 
the B2 background galaxy has My ~ -22.0 and the blue colors 
typical of a star-forming galaxy. 




Fig. 6. -TheB2/B3 system of arcs (see text for more details). Cluster eUipticals #46, #76 
and #23 are also labelled. Each tick mark represents V , and the orientation is as in Fig. 1 . 
The dotted lines indicate the path of the z = 3.5 critical hne close to the observed arcs. 

B3 consists of three images; B3b/c is adjacent to cluster el- 
liptical #23, and B3a, the counter image, lies outside the critical 
line, ^ 1" below and to the right of B2c in Fig. 6. B3d appears 
to be a singly-imaged component of the background galaxy that 
produces B3a/b/c. This interpretation is supported by the pho- 
tometric analysis in Table 1. We also estimate that B3 is a star 
forming galaxy with a similar luminosity (My ^ -22.4) to B2. 

B4 (Fig. 1) — This is another complex and unusual arc, the 
precise nature of which is uncertain. One possibihty is that 
B4 is part of a multiple image system containing B6; however, 
photometric analysis (Table 1) does not support this hypothe- 
sis. More complete color information will improve our ability 
to comment on this arc. If B4 is a singly imaged source it lies 
at a redshift of z ^ 1 -0. 

3.1.3. Faint Arcs 

B5 (Fig. 1) — B5 is a very faint arc lying in the saddle re- 
gion between the central galaxy and cluster elliptical #99. We 
interpret this arc as an image of a background galaxy which is 
magnified by this saddle and therefore constrain the redshift of 
the galaxy to be z ^ 2.4. 
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B9a & B9b (Fig. 1) — This is a pair of extremely faint arcs 
detected only because of the superb resolution and low sky 
background of the HST observations. Both arcs appear to com- 
prise three images resulting from the magnifying effect of the 
saddle potential between the cluster center and cluster elliptical 
#260. We estimate that they lie at a redshift of z > 3. 

B14 (Fig. 1) — B14 is barely detected above the sky back- 
ground in the HST frame, however we detect it at A" = 19.70 ± 
0.05 in the UKIRT frame, giving (R-K) = 6.02 ±0.11. As 
B14 is singly-imaged, we obtain a redshift limit of z < 3.9 from 
our lens model. A full analysis of galaxies exhibiting extreme 
colors and the general background population lying behind our 
cluster sample will be presented in a future paper. 

3.1.4. Arclets 

We also detect numerous singly-imaged arclets in the field of 
A 383 (Fig. 1 and Table 1) and have obtained upper limits on 
the redshifts of three of these arclets (B6, B7 and B8) on the 
basis of their singly-imaged nature. 

B6 (Fig. 1) — We obtain a redshift limit of z < 0.9 for this 
arc let, however it's blue (B-R) and (R-I) colors (Table 1) are 
consistent with it being a star forming galaxy at z ~ 1 .5-2.0. 
The latter interpretation supports the possibility that B6 and 
B4 contain multiply imaged features of a single background 
galaxy. As noted in §3.1.2, improved photometry of B4 will 
help to resolve this uncertainty. 

B7 (Fig. 1) — B7 lies at a redshift of z < 1.3 which, on the 
basis of its (B-R) color (Table 1) is consistent with B7 being a 
moderately star forming galaxy. 

B8 (Fig. 1) — This arclet has a redshift limit ofzS^ 2.0, while 
its colors are broadly consistent with this being a star forming 
galaxy at z ~ 1.0-1.4. 

3.2. Lens Model 

The mass distribution of A 383 was reconstructed using the 
parametric lens inversion method described in detail by K96. In 
addition to a cluster-size mass component, the model contains 
29 individual galaxy masses (including the central galaxy) cov- 
ering the entire HST field of view. The parameters describing 
the model were constrained were constrained within the cluster 
core using four of the multiple image systems identified in the 
HST frame (BO, Bl, B2 and B3) and on larger scales using the 
weak shear field estimated from the shapes of the faint galax- 
ies across the whole frame. The cluster scale mass that dom- 
inates the mass distribution has a core radius of (46 ± 3) kpc, 
velocity dispersion of (900 ± 20)kms"' and a cut-off radius of 
(1200 ± 100) kpc (see K96 for details of how these quantities 
are defined). 

The mass distribution obtained from the lens model is over- 
laid on the optical image in Fig. 1 . We calculate the total pro- 
jected mass within the radius of the giant arc (65 kpc) to be 
(3.5 ±0.1) X IO'^Mq. Using the weak shear measurements we 
can extend our model to larger scales and constrain the pro- 
jected mass within a 250 kpc radius of the cluster center to be: 
(1.8 ±0.2) X IO'^Mq. 

In Fig. 7 we plot the projected surface mass density derived 
from the lens model as a function of the distance from the 
cluster center (r). On scales r > 10 kpc, the density profile is 
dominated by the cluster-scale potential, the core radius of the 
cluster-scale mass distribution 46 kpc) being consistent with 
previous estimates of core radii from lensing studies (Mellier et 
al. 1993; Small et al. 1996; Natarajan et al. 1998). The shape of 
A383's profile is however only constrained out to r ~ 400 kpc. 



this being the extent of the WFPC2 field of view. The mass 
profile is most accurately measured within r ^ 100 kpc where 
information from the multiple images is available. 



Radius (kpc) 

1 10 100 




1 10 1 00 



Radius (arcsec) 

Fig. 7. - The projected surface mass density (S) of A 383 as a function of radius, obtained 
from the lens model by measuring the mass density in concentric annuli. The dashed lines 
illustrate the difference between the profile slope at the positions of the two radial arcs 
(BOb and Bid), which are marked by arrows at 6.0 and 19.9 kpc respectively. 



3.2.1. Radial Arc Properties 

BOb and Bid (see Figs. 1 & 5) lie at distances of (6.0 ± 
0.4) kpc and (19.9 ± 0.4) kpc from the cluster center respec- 
tively. These two arcs therefore enable us to probe the shape 
of the mass distribution on scales of r < 20 kpc. Specifically, 
the radial arcs constrain the slope of the density profile, as their 
radial positions depend on the local gradient of the projected 
mass, rather than on the mean enclosed surface density. 

A detailed analysis of the shape of A383's density profile is 
currently limited by the pseudo-isothermal elliptical functional 
form adopted in the lens model (see K96). However, in or- 
der to estimate the slope of the density profile at the positions 
of the two radial arcs, we fitted various polynomial functions 
to the density profile data. De-projecting this analysis to ob- 
tain the slope of the three-dimensional density profile (p), we 
obtain the gradients: dQog p) / dQog r) = -1.41 and -1.29 at 
r = 6.0kpc (BOb) and r = 19.9 kpc (Bid) respectively; the error 
on these slopes is of the order < 0.03. 

We interpret this difference in slope at the positions of the 
two radial arcs, specifically the steeper slope at r ^ 6 kpc, as be- 
ing due to the presence of a distinct central galaxy halo within 
the overall cluster halo. 

Williams et al. (1999) suggest that radial arcs are a useful test 
of the lensing role of such central galaxy halos within an over- 
all cluster-scale NFW halo. They predict that \og(0r/ 9t) (where 
9r and 9, are the angular positions of the radial and tangential 
arcs respectively) should decrease with decreasing cluster mass 
as the lensing role of the central galaxy becomes more signifi- 
cant relative to the overall cluster mass distribution. From their 
Fig. 7 we estimate that for A 383: \og{9r/9,) = -0.35 to -0.65, 
taking into account the velocity dispersions obtained from spec- 
troscopy (§2.3) and the lens model (§3.2). 
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Table 1 

Pholomelry and Speclroscopy of Lensed Fealures in A 383 



Feature 


R 


(B-R) 


(R-I) 


a-io 


z 


nd 

source 


BOa 
BOb 


21.98±0.01 
21.95 ±0.24 


1.02 ±0.01 


0.74 ±0.01 


1.91 ±0.03 


1.0103 ± 0.0001° 


25.0 ±0.3 


Bla/b 
Blc 

Bid 


21.50±0.02 
23.94 ±0.07 


1.14 ±0.03 
1. 1 1 ± 0.03 


1.07 ±0.04 
1 no ± 03 


2.06 ±0.14 
1 qn ± 07 


1.1±0.1 


26.7 ±0.6 


B2a/b 

Blc 

B2d/e 


23.79 ±0.06 
23.90 ±0.01 
22.56 ±0.21 


0.87 ±0.03 
0.75 ±0.26 


0.29 ±0.02 
0.54 ±0.23 


1.94 ±0.08 
< 2.3 


3.5 ±0.2 


26.8 ±0.7 


B3a 

B3b/c 

B3d 


23.80 ±0.02 

ZZ.oZ ± U.Z4 

23.86 ±0.02 


0.83 ±0.02 

U.6D ± 0.2/ 

1.09 ±0.03 


0.40 ±0.02 

A <0 -I- A 0"7 

0.50 ±0.02 


< 1.9 

< 2.5 

< 2.0 


3.5 ±0.2 


26.3 ±0.2 


B4 


22.07 ± 0.01 


0.94 ±0.26 


< 0.4 




S, 1.0* 




B5 


23.35 ±0.03 








> 2.4 




B6 


23.42 ±0.01 


0.42 ± 0.03 


0.12 ±0.07 


< 3.5 


S,0.9'' 




B7 


24.23 ±0.02 


0.63 ±0.02 


<0.5 




S, 1.3* 




B8 


24.68 ± 0.05 


0.84 ±0.14 


0.81 ±0.17 


< 4.1 


& 2.0* 




B9a 
B9b 


25.20 ±0.11 
24.48 ±0.08 








> 3.0 

> 3.0 




BIO 


24.85 ±0.05 


0.19 ±0.03 


0.57 ±0.07 


<4.5 






Bll 


22.92 ±0.01 


1.60 ±0.01 


1.02 ±0.01 


2.46 ±0.01 






B12 


23.56 ±0.02 












B13 


23.97 ±0.01 


1.45 ±0.14 


1.05 ±0.16 


2.19±0.18 






B14' 


26. 12 ±0.02 


> 1.2 


<2.0 


>4.0 


S,3.9 




B15 


22.23 ±0.02 


1.35 ±0.01 


0.94 ±0.02 


1.96 ±0.09 


1.1±0.1 




L, Cluster E/SO 


18.01 


2.44 


0.67 


2.50 


0.188 





a) Spectroscopic redshift; b) Redshift assumes single image; c) Extremely Red Object with (S - A") = 6.02 ± 0. 1 1 . d) Estimated iJ-band magnitude of the background galaxy, corrected for 
lens amplification; the error bars reflect the uncertainties in the lens model and generally scale with the number of multiple images per backgrotmd galaxy. 



We have computed angular position ratios of log{6r/6,) = 
-1 .04 ±0.07 and -0.5 1 ±0.03 for A 383's inner (BOb) and outer 
(Bid) radial arcs respectively. The outer radial arc appears to be 
consistent with Williams et al.'s predictions, however the inner 
radial arc falls well outside their prediction. This discrepancy 
may be due to simplifying assumptions made by Williams et al. 
regarding the nature and complexity of cluster substructure. We 
also note that, due to the (almost) circular symmetry of A 383, 
the radial arcs are not forced to lie at the radial critical line in 
the image plane. This is another potential source of discrepancy 
between the observations and the predictions. 

3.2.2. Comparison With Numerical Simulations 

Ghigna et al. (2000) have performed numerical simulations 
of galaxy clusters at the highest resolution to date, resolv- 
ing structure down to scales of r ~ 5 kpc. They find a de- 
projected slope of -1.6 ±0.1 on scales of r = 5-100 kpc in a 
(J ^ 800 km s"' cluster This is steeper than the density profile 
of A383 on comparable scales, ~ -1.3 (see §3.2.1 and Fig. 7). 
We will present a detailed comparison of the lens model using 
NFW and Moore et al. (1998) density profiles in a subsequent 
paper. 

Recent progress has also been made in the inclusion of gas in 
numerical simulations (e.g. Pearce et al. 1999). However, the 
resolution of these simulations precludes analysis of the form 
of the mass profile on the scales probed by the radial arcs in 
A 383 (r < 20 pc). Higher resolution simulations including gas 
cooUng have been performed by Lewis et al. (1999). However, 
these authors concede that the mass of the central galaxies in 
their simulations is unrealistically high, which hampers a reli- 



able comparison with our observations of A 383. 

The constraints on the density profile derived from the ra- 
dial arcs in A 383 confirm the need for a massive central galaxy 
to explain the lensing properties of this X-ray luminous clus- 
ter. Our observations therefore support the case for increasing 
the resolution of N-body simulations of clusters including gas 
physics. 

3.3. X-ray Properties of A 383 

The radial X-ray surface brightness profile around the cluster 
X-ray centroid is shown in Fig. 8. We attempt to describe the X- 
ray emission profile by fitting a standard beta model (CavaUere 
& Fusco-Femiano 1978) 

/w=/o(l+(4)Y'^^■^ 

convolved with the HRI point-spread function (PSF), to the ob- 
served X-ray surface brightness distribution. The free fit pa- 
rameters are the peak surface brightness (/o), the 3-d core ra- 
dius ('"c''), and the slope parameter (/3). We note that the use 
of the beta model implies that the gas is isothermal and in hy- 
drostatic equilibrium. While these assumptions are likely to be 
violated in most clusters, numerical simulations show that us- 
ing the beta model introduces, on average, no bias and only a 
moderate scatter of typically 20% in the cluster mass estimates 
(Evrard, Metzler & Navarro 1996). 

We fit two beta models, one to the full radial range of the 
data including the likely cooling flow region (model A), and a 
second one to the outer regions only (model B). While model 
A yields a poor description of the X-ray emission at large radii 
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(see Fig. 8) it provides us with an acceptable fit to, and an an- 
alytic description of, the emission profile at r < lOOkpc, which 
can be used to derive the gas density and, assuming hydro- 
static equilibrium, the total binding mass within the core re- 
gion (see below). Model B, obtained by fitting only the range 
from r= 12kpc to 850kpc (3.5'), describes the non-cooling- 
flow component (Fig. 8), for which we find r^'' = (23 ± 3)", 
corresponding to (93 ± 12)kpc, and (3 = 0.65 ±0.03, consistent 
with the canonical value of 2/3. The quoted errors include sys- 
tematic uncertainties which we explored by varying the radial 
range used for the fit. 



radius (kpc) 




10 

radius (arcsec) 

Fig. 8. - The radial X-ray surface brightness profile of A 383 as observed with the ROSAT 
HRI. Overlaid are the best-fit beta models (convolved with the HRl PSF) and fitted to the 
full profile (model A, dashed line) or to the outer regions only (model B, solid line). The 
excess emission in the cluster core is caused by an additional cooling flow component, 
which is not included in the models. 

Subtracting beta model B from the observed profile we find 
the cooling flow to contribute 13% to the total cluster flux 
within a radius of 500 kpc, corresponding to a luminosity of 
1 X 10** ergs"' (0. 1-2.4 keV). A rough (usually low) estimate 
of the mass deposition rate, M, can be obtained by converting 
this excess luminosity to the bolometric band (by multiplying 
by a factor of 3.03 appropriate for the assumed cluster gas tem- 
perature of 7.1 keV) and equating it to the luminosity of the 
cooUng flow 

T ^ UT 

^cool — _ 

2 /iWp 

(Fabian 1994). We thus derive a lower limit to the mass de- 
position rate of M > 200 yr"\ a value typical of moderate 
cooUng-flow clusters (e.g. Peres et al. 1998). 

Under the assumption of hydrostatic equilibrium, and using 
a beta model to describe the radial gas density profile, the total 
gravitational mass within a projected radius r can be derived 
from 



M« r) : 



AtmpG 2 ^l + (r/r3d)2 



(Soucail et al. 2000). 

Using the parameters of model A we estimate the projected 
mass within a radius of 65 kpc to be 4.0^{ j x \0^^;Mq, in good 
agreement with our lensing mass measurement of (3 .5 ± 0. 1) x 
lO'^M© within the same projected radius (§3.2). The equiv- 
alent mass estimate within a radius of 250 kpc is more uncer- 
tain, but is approximately (1.2 ±0.5) x XO'^'^Mq, again in rea- 
sonable agreement with the lensing measurement. We note that 
the bulk of the uncertainty in the X-ray mass estimates is caused 
by the adopted error of ±2 keV in the global gas temperature. 



kT . Forthcoming observations of our full cluster sample with 
Newton will allow us to obtain much more accurate X-ray mass 
measurements on scales of ~ 10- 1000 kpc. 

4. DISCUSSION AND CONCLUSIONS 

We have discovered numerous new gravitationally lensed 
features, including a giant arc and two radial arcs in A 383, 
a massive, X-ray luminous cluster of galaxies at a redshift of 
z = 0. 188. This is the first cluster in which two radial arcs have 
been identified and these provide a detailed view of the cluster 
mass distribution on scales of r < 20 kpc. 

The morphologies and positions of the arcs have been used, 
in conjunction with ground-based multi-color photometry, to 
constrain the mass distribution and compute a precise mass den- 
sity profile of this cluster within a radius of ~ 250 kpc. 

Within the core radius (r = 46 ± 3 kpc) derived from the 
lensing analysis, the density profile exhibits a shallower slope 
{d{\ogp)/d{\ogr) = -1 .29 ± 0.03) than a Moore profile (Moore 
et al. 1998; Ghigna et al. 2000) and a steeper slope than a single 
NFW profile. 

We have also analysed the properties of the two radial arcs, 
interpreting the difference in angular position of these arcs as a 
signature of the lensing role of the massive central galaxy. This 
supports the proposal by Williams et al. (1999) that massive 
central galaxy halos are required to explain the lensing proper- 
ties of the cores of massive clusters (a ^ 1000 km s"'). How- 
ever, only the outer radial arc of A 383 falls within the range 
predicted by Williams et al. for the the angular position ratio of 
radial and tangential arcs. The inner radial arc lies well outside 
of their predicted range. This discrepancy is probably caused by 
simplifying assumptions in Williams et al.'s theoretical models. 

The lensing analysis is complemented by an analysis of the 
cluster's X-ray properties as obtained from archival ROSAT HRI 
data. We find strong evidence for a cooling flow in agreement 
with the optical properties of the cluster, and derive a lower 
limit on the mass deposition rate of 200 Mq yr"' , typical of 
moderate cooling flow clusters. The X-ray estimates for the to- 
tal mass within the projected radius of the giant arc are in good 
agreement with the lensing measurements and are consistent 
with earher claims (e.g. Allen 1998) that significant discrepan- 
cies between lensing and X-ray mass estimates are caused by 
substructure in unrelaxed, non-cooUng-flow systems. 

We note that the cluster sample used in our survey is effec- 
tively X-ray selected, thus allowing us, for the first time, to 
measure in detail the mass distribution in massive clusters of 
galaxies unaffected by selection biases. The theme of our fu- 
ture program will be to analyse the lensing signal and X-ray 
emission (as observed with Newton) from all twelve clusters in 
our sample and to trace the form of the mass profile from 50 kpc 
to 5 Mpc. We wiU investigate the dispersion in the profiles and 
the core properties of the clusters and determine how these cor- 
relate with the cluster's dynamical state and the presence of a 
cooling flow. We shall also attempt to constrain the high end 
of the cluster mass function and calibrate directly the cluster 
mass-temperature and mass-luminosity relations at z ~ 0.2. 
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